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We use data from the WMAP temperature maps to constrain a scale-dependent generalization 
of the popular 'local' model for primordial non-Gaussianity. In the model where the parameter 
/nl is allowed to run with scale k, /nl(&) = /NL(' c /fcpiv) n ' fNL , we constrain the running to be 
n /nl = 0.30*i'§ at 95% confidence, marginalized over the amplitude f^ L . The constraints depend 
somewhat on the prior probabilities assigned to the two parameters. In the near future, constraints 
from a combination of Planck and large-scale structure surveys are expected to improve this limit 
by about an order of magnitude and usefully constrain classes of inflationary models. 
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Introduction. Non-Gaussianity in the distribution of 
primordial density fluctuations provides a unique window 
into the physics of inflation. The magnitude of primordial 
non-Gaussianity and its dependence on scale provide in- 
formation about the dynamics of scalar field(s), their in- 
teractions, and the speed of sound during inflation. Con- 
straints on non-Gaussianity have traditionally come from 
the measurements of the three-point correlation function 
of the cosmic microwave background (CMB) temperature 
anisotropies. Upper limits from COBE pQ have been im- 
proved by two orders of magnitude by the WMAP ex- 
periment [2] . Moreover, clustering of galaxies and galaxy 
clusters has also been identified as a powerful probe of 
non-Gaussianity [5], already leading to interesting con- 
straints that are complementary in their information con- 
tent to the CMB measurements. 

So far most attention has been devoted to the "lo- 
cal" model of primordial non-Gaussianity, where the pri- 
mordial Newtonian potential 4>(x) is modified with a 
quadratic term: 4> — 4>G + /nl(</>g — ( < Pg))i wnere 4>G is 
a Gaussian potential [1] . The parameter /nl is currently 
constrained to be 32 ±21 by WMAP (0; see also [51 1]) 
and 28 ± 23 by the large-scale structure [7HS]- Several 
other non-Gaussian models have been constrained as well 
(e.g. pTTlfTTj ). However, the 'running' with physical scale 
of these models, which may carry important information 
about the number of inflationary fields and their interac- 
tions [12 22], has not yet been constrained with current 
data (except for a very rough estimate of the angular- 
multipole dependence of /nl [H] and implicit constraints 
on a braneworld-motivated model 23 ) . Such constraints 
have only been forecasted for future experiments |24H28j . 
Constraining the running of non-Gaussianity therefore 
presents a major new opportunity to probe inflationary 
physics, and is just becoming feasible. In this Letter, we 
present the first such constraints. 

Model. In this work we consider a physically motivated 
generalization of the local model, where the parameter 
/nl is promoted to a function of scale k. In particular, we 
seek to constrain the two-parameter power-law subclass 



of the generalized models [25] 



/nlO) = /nl 



Vnl 



(1) 



where k* is an arbitrary fixed parameter, leaving / NL 
and n/ NL as the parameters of interest in this model. 
Such scaling is expected in inflation when more than one 
field dominates or when there is self-interaction, and its 
signatures in the CMB and LSS have been discussed in 
the literature j24j [25j [29] . The parameter n / NL is often, 
though certainly not always, expected to be < 0(1) in 
inflationary models, but in our phenomenological model 
it is allowed to take any value. 

Bispectrum and / N l estimator. The primordial bispec- 
trum of the fNh(k) model from Eq. ([!]) is straightforward 
to calculate: 

F(ki, k 2 , h) = 2 [f N MP(k 2 )P(k 3 ) + perm.] , (2) 

where the full bispectrum is B(k\,k2,kz) = (27r) 3 (5(fci + 
k2+k^)F{k\, ki, k%). Here P is the power spectrum of the 
primordial curvature perturbations, and 5 is the Dirac 
delta function. 

Constraining the running parameter n/ NL seems diffi- 
cult because of the apparent requirement to find an esti- 
mator for a parameter in an exponent. To avoid this, we 
resort to an indirect approach where, for a fixed value of 
n / NL , we estimate the parameter / NL using modifications 
of the well-known KSW estimator [30], which is known 
to be nearly optimal [3U 131] • We then iterate over the 
values of the running ?i/ NL to obtain the full likelihood 

^(/nL^/nl)- 

The theoretical expectation for the bispectrum of the 
temperature anisotropies in the cosmic microwave back- 
ground can be explicitly evaluated, starting from the def- 
inition of the generalized non-Gaussian local model in 
Eq. ([!]) to account for the running n/ NL : 

□theory//.* \ r> r* r ., 

UND^/nlJ = ^NL«i<2^ x 
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r 2 dr (a £l (n/ NL , r)ft 2 {r)(i h (r) + perm.) (3) 



2 



where is the Gaunt integral and 

2 1 



ott(r) 



" K piv 



fc 2 +»/ NL t e (k)j £ {kr)dk (4) 



ft(r) = ~ / eP^k)U{k) H (kr)dk. 

7T ' 



(5) 



Here, tt is the radiation transfer function, which can be 
calculated using CAMB [33]. Following KSW [30] we can 
define new, filtered maps A(h, r) and B(ri, r), 

A(h,r) = at (n/ NL , r ) ^ Of mXtm (&) , (6) 

B(n, r) = y^j8f(r)^-Of m yg m (fi). (7) 

Then, we write down the skewness 5(ra/ NL ): 

S(n fN J = J r 2 dr J d 2 hA(n,r)B 2 (h,r), (8) 

which requires n/ NL as input (through A), and does not 
require a priori knowledge of / NL . 



The observed CMB bispectrum is defined as B^f' e 

therefore reduces to 



), and S(ra/ NL 

obs Stheory 



E 



^'xtits^iitita ^ NL _ ^ 



C^ 2 Q. 



(9) 



theory 



and 6/) is the beam 



where = b^b^b^B , 

transfer function. 

We now define F = F(n/ NL ), the Fisher matrix for 
/nlj equivalent to the cumulative signal-to- noise squared 

1 



of the theoretical bispectrum for / NL 

/ Athcory / r 



n^ NL )= E 

lx<lz<l 3 



NL 



1) 



C^i w 2 Cf 3 



The theoretical expectation for 
cubic KSW estimator for / NL is: 



/nl — 



S_ 
F' 



oc / NL , so the 



(11) 



We used HEALPix, by way of HealPy, to do the for- 
wards and backwards spherical harmonic transforms re- 
quired to obtain the A and B maps. 



Cut-sky maps. Equation (111 works well for a full-sky 



map, but a sky cut introduces a spurious non-Gaussian 
signal. To account for the masking of the CMB sky, we 
make the substitution S — > S cu t = S/ / s k y + Senear [HI] . 
Siincar is an addition to skewness from Eq. Q, calibrated 
to account for partial-sky observations: 



1 

/sky 



2 dr J d 2 n[A(n ir )(B 2 im (n,r)} 



MC 



+25(n, r) (A sim (h, r)£ sim (n, r)) 



MC\ 



(12) 



The subscripted filtered maps A s 
from Python-produced Monte Carlo realizations of the 
cut CMB sky; the brackets Que indicate an average over 
all 300 Monte-Python maps. The simulated maps were 
produced using the prescription laid out in Appendix 
A of the WMAP5 paper [35]; the only difference (aside 
from using the WMAP7 cosmological model) is that we 
used a uniform weighting for the maps, rather than the 
slightly more complicated weighting given there, since it 
only gives a marginal improvement in estimating /nl- 

Likelihood Evaluation. To find the likelihood, we first 
find a x 2 statistic for / NL , given a value of nt Nli . Taking 
the angular-averaged bispectrum Be 1 g 2 i 3 as our observ- 
ables, we have: 

X 2 (/nL^/nl) = 

^ (BfMata ~ /nl^iI^^/nd/nL = 1)) 

(Again, this works because the theoretical expectation for 
Bi 1 £ 2 i 3 oc /nl-) Using Eqs. Q and (10), we can rewrite 
X 2 as 



x 2 (/, 



NL) w /nl 



,)=F / NL - 



xl 



§1 
F ' 



(14) 



the 



where xl = E Ws ( B tU) I^C l? C lz ) is 
goodncss-of-fit parameter for the data with respect to 
the /nl — case. Note that the numerator of Xo is an 
observed quantity, and the denominator is based solely 
on the theoretical prediction for the power spectrum (as 
well as a few noise and beam parameters of WMAP). 
Therefore, Xo does not depend on / NL or n / NL at all. We 



can use the definition of / NL in Eq. (11) to rewrite the 
expression for x 2 a s follows 

( 10 ) X 2 (/nl,«/ n J=^(/nl-/nl) 2 + x^-(/nl) 2 ^ (15) 



the x 2 

— JNL' 



is, as expected, 
and one obtains 



For a fixed value of n / NL , 
minimized in / NL when / NL 

A more interesting task is to calculate the constraints 
when ri/ NL is allowed to vary. With an expression for 
X 2 (/NL' n /NL) m hand, we can write an expression for 



the likelihood, £(/ NL , nj 
constant term with Xo) 



ex exp(— x /2) (dropping the 



£(%l./nl) oc exp 



F[f*(L-fi 



'NL 



exp 



(f^) 2 F 



(16) 



To marginalize over / NL is also straightforward 



c < " Jsl 1 = / ^KnlJnl) dfixh oc -^ exp 



(f^) 2 F 



(17) 
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FIG. 1. Likelihood in the ?i/ nl -/nl plane (left panel) and marginalized over (right panel). The principal constraints, 
shown in the left panel and with the bold blue curve on the right, correspond to the flat prior on /^l at the pivot value where 
the constraints on /n L and fi/ NL are uncorrelated (see Eq. (19l). In the right panel we also show the marginalized likelihood 
for n/ NL with a prior on f^ L that is uniform in log(/]JJ L ) f° r |7nlI > 0.1 and zero otherwise. The dashed curve in the left panel 
shows the quantity /nl> which is the best-fit value of the parameter for a fixed «-/ NL - See text for other details. 



''/n 



and k. 



where, recall, -F(?i/ NL ) is denned in Eq. ( |10[ ). 

WMAP7 constraints on n/ NL . Figure [l] shows the like- 
lihood C in the n/ NL - plane, as well as the likelihood 
for n/ NL alone, calculated from the WMAP7 temperature 
maps. We used a weighted and masked combination of 
the WMAP V and W band maps with the monopole and 
dipole subtracted, as recommended by the WMAP team 
[55] . To extract full information from WMAP maps, 
we used multipoles out to £ max = 800 for the sums in 
Eqs. (Rjj), ([7]) and (10). We did not find a significant im- 
provement between £ max — 700 and £ max — 800; we chose 
the higher value to be conservative in our analysis. 

The quantity x 2 is independent of our choice for fc p i v , 
but the likelihood itself is not, since F is inversely pro- 

2n /■ 

portional to k iv . The true pivot scale favored by the 
data is the value of k p i v for which the errors in are 

uncorrelated with the errors in n y NL . We find this scale a flat prior on /j^ L and 



shown at the left panel of Figure jl] assuming flat priors 
on and 



= k^ Atr} ps 0.064 /iMpc" 



Dependence on the prior. As with most present-day 
cosmological measurements, the precise constraints de- 
pend on the prior probability on the parameters we are 
constraining. Even for a simple flat prior on and 
rt/ NL , the actual effective prior depends on the a priori 
chosen pivot in wavenumber fc*. For example, a flat prior 



on (/i 



NU 



1(1) 



by using the likelihood to calculate the covariance matrix 
C between and n< 



/NL(fc*,i) defined at some pivot corre- 
sponds to a non-flat prior on some (/nl)' 2 '' = /nl(^*,2) 
defined at some other pivot £^.2, since (/nl) 1 ' 2 ' = 
(/nl)^ 1 H^*,2/&*,i)" /nl • If we assume some alternate 
pivot 2 but hold the flat prior in /^_r , the contours 
in the n/ NL -/^ L plane (left panel of Fig. IT]) are stretched 
vertically by a factor of (£^2/0.064 h Mpc" 1 )"^ . 

We have experimented with different k-pivot values for 
i/ NL . We have also investigated 



C p*,p 3 = ((Pi-Pi)(Pj ~Pj))- 



(18) 



We can easily find the pivot value fc p i v that diagonalizes 
the covariance matrix C (see e.g. Ref. [26] ) 



/c* exp 



other possibilities, such as the prior that assigns equal 
weight to each decade in above 0.1 (so uniform in 

log(f^ L ), but cut off at the arguably lowest-ever observ- 
able value of I/nlI = 0-1 so that the total integrated 
likelihood is finite). We present the two aforementioned 
examples, showing constraints on fi/ NL marginalized over 
/^lj m the right panel of Fig.[TJ In the end, we decide to 
quote results for the flat prior and the uncorrelated fc p j v 



fm^n 



(19) 



value from Eq. (19), which most closely follows priors to 



where fc* is the (arbitrary) pivot used initially, and 
is the corresponding value used in C. Despite the fact 
that fc* and show up in the expression, k plv does 
not depend on them: it is a fixed number telling us 
roughly where the experiment has greatest power (and 
where normalization and running of fisihik) are precisely 
uncorrelated). We find that k^ AP7 ps 0.064 h Mpc" 1 . 
The 68%, 95%, and 99% constraints on /j^ L and n/ NL are 



both non-Gaussian and other cosmological parameters 
applied in the literature. 

Putting it all together, we can get the estimate for 
ri / NL from the WMAP7 data for a flat prior on at 
the pivot /c piv from Eq. The 68% (95%) confidence 

interval is 



0.30 



+0.78 (1.9) 
-0.61 (1.2) 



(20) 



The current constraints are therefore fully consistent 
with no running, as Fig. [TJ clearly indicates. Figure [2] 
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FIG. 2. Constraints propagated to /nl(&)- We also show 
several models that are reasonable fits to the data (all within 
the 99% confidence limit of the left panel of Fig. [T]) to guide 
the eye as to how typical models from our ansatz behave. 



shows the constraints in the /nl(&) plane together with 
a few representative models allowed by the data. 

Conclusions. We have presented the first constraints 
on the scale-dependence of (any form of) non-Gaussianity 
using the WMAP7 data. The constraints are compatible 
with zero running, ny NL = 0, with very mild (< 1-sigma) 
preference for a positive value of n/ NL . We will learn 
more soon: the Planck data and the data from upcoming 
large-scale structure surveys should be able to improve 
constraints on the running of non-Gaussianity by about 
an order of magnitude (2H [37] HI] , thus shedding impor- 
tant new light on the physics of inflation. 
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